The spin Hall effect as a probe of nonlinear spin fluctuations 
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The spin Hall effect and its inverse play key roles in spintronic devices since they allow conversion 
of charge currents to and from spin currents. The conversion efficiency strongly depends on material 
details, such as the electronic band structure and the nature of impurities. Here we show an anomaly 
in the inverse spin Hall effect in weak ferromagnetic NiPd alloys near their Curie temperatures with 
a shape independent of material details, such as Ni concentrations. By extending Kondo's model 
for the anomalous Hall effect, we explain the observed anomaly as originating from the second-order 
nonlinear spin fluctuation of Ni moments. This brings to light an essential symmetry difference 
between the spin Hall effect and the anomalous Hall effect which reflects the first order nonlinear 
fluctuations of local moments. Our finding opens up a new application of the spin Hall effect, by 
which a minuscule magnetic moment can be detected. 



PACS numbers: 

Spintronics has provided an impetus for the intensive 
development of techniques to generate or detect spin cur- 
rents, such as nonlocal spin injection^, spin pumping- 2 
and spin light emission devices 3 -. The spin Hall effect 
(SHE) 4 , which converts a charge current to a spin cur- 
rent, and its inverse have been widely explored in semi- 
conductors 3 -"— and nonmagnetic metals 2 -^ by using such 
techniques. While most of the work has concentrated 
on the phenomenon and mechanism of the SHE itself, 
Uchida et al. have used the SHE as a detector of a spin 
dependent chemical potential arising from the spin See- 
beck effect^. Further applications of the SHE are, how- 
ever, still lacking. Recently the doping of gold with in- 
dividual iron impurities was considered as a mechanism 
to enhance the SHE through spin fluctuations of the iron 
moments 1 ^. This raises the question whether the SHE 
might be sensitive to collective fluctuations, in particular 
close to a magnetic phase transition where the fluctua- 
tions are strong. As for the anomalous Hall effect (AHE) 
in a pure ferromagnetic metal such as Fe 1 - 1 - and N i 11 ' 12 , a 
sharp peak of the Hall resistivity (ph) appears only below 
the Curie temperature (lb). This can be understood as 
an effect of fluctuations of local magnetic moments 1 ^. In 
this report we shall demonstrate experimentally, and pro- 
vide a theoretical argument, that the inverse SHE (ISHE) 
also shows an anomaly, but that it is substantially dif- 
ferent from that of the AHE. Figure 1 shows schematic 
images of the AHE and the ISHE of a weak ferromag- 
netic material near lb. In the AHE, a charge current 
is injected into the ferromagnetic material, and the dif- 
ference of the skew scattering probabilities for spin-up 
and spin-down electrons results in a Hall voltage. In the 
ISHE, on the other hand, a spin current is injected and 
both spin-up and spin-down electrons are deflected to 



the same side, which also results in a Hall voltage. Con- 
trary to the AHE, however, we observe a dip and peak in 
the ISHE for weak ferromagnetic NiPd alloys below and 
above lb, respectively. Since the quantity of the involved 
Ni ions are very small, the present system could be used 
for sensitive magnetometry. 

Results 

AHE and ISHE of NiPd alloys. 

We measured the SHE of a weakly ferromagnetic NiPd 
nanowire with 7%, 8%, and 9% Ni using the nonlocal 
spin injection technique^ 4 -"— at 5 — 50 K. Figures 2a 
and 2b show the schematic structure and scanning elec- 
tron microscopy (SEM) image of the lateral spin valve 
with the NiPd middle wire. The NiPd alloy is a weak 
ferromagnet when the Ni concentration is larger than 
2.3%i£. Depending on the Ni concentration, its Tq varies 
from a few K to hundreds of K. In order to charac- 
terize the magnetic property of NiPd wires in our de- 
vices, AHE measurements were performed for the refer- 
ence NiPd Hall bars prepared at the same time as the 
SHE devices. A representative AHE resistance R xy for 
Nio.08Pdo.92, which is proportional to the magnetization, 
is shown in Fig. 3a. We note that in the AHE measure- 
ments, the external field H was applied perpendicularly 
to the plane, which is different from the ISHE measure- 
ments (see Fig. 2b). At T = 5 K a hysteresis loop of the 
magnetization can be clearly seen. With increasing tem- 
perature, the magnetization continually decreases to zero 
and R xy simply shows the ordinary Hall resistance above 
30 K. The Curie temperature lb for Nio.osPdo.92 can be 
determined to be about 21 K by plotting the remnant 
magnetization M r and coercive field H c versus tempera- 
ture (Fig. 3b) , which is lower than the value for the bulk 
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FIG. 1: Principles of the anomalous Hall and inverse 
spin Hall effects near Tc in a weakly ferromagnetic 
metal, a, AHE (incident charge current Jc) and b, ISHE 
(incident spin current Js). The green and pink spheres rep- 
resent the conduction electrons (i.e., — |e|) and localized ones, 
and the grey arrows show the directions of spins. The fluctu- 
ation of the localized spins near Tc is indicated by the longer 
grey arrows with shades. The skew scattering probabilities 
for spin-up and spin-down conduction electrons, resulting in 
anomalous Hall or inverse spin Hall voltage, are represented 
by red and blue arrows respectively. 



alloj*i£. This is presumably due to the imperfect distri- 
bution of Ni in our alloy, prepared by the ion implan- 
tation technique, compared with bulk alloys (see Meth- 
ods). We also performed the same measurements for the 
other two concentrations where Tc for Nio.07Pdo.93 and 
Nio.09Pdo.91 arc 16 and 32 K, respectively. 

The ISHE is based on the spin to charge current 
conversion, which follows the vector product as below: 
Jc oc Js x sc where Jc and Js are the charge and spin 
current density, sc is the spin orientation of the con- 
duction electron. Our device configuration is depicted 
in Fig. 2a. The charge current flows from one of two 
Nio.8iFeo.19 (hereafter Py) wires, i.e., Pyl to the left hand 
side of the Cu wire. Since the spin-orbit (SO) interac- 
tion of Cu is quite weak, it works as a transmitter of 
the generated spin current at the Cu/Pyl interface to 
the right hand side of the Cu wire. The spin current is 
preferentially absorbed into the NiPd middle wire below 
the Cu wire because of its strong SO interaction. This 
absorption is confirmed by a reduction of the nonlocal 
spin valve (NLSV) signal detected at the second Py elec- 
trode Py2, as detailed in refs. 15 and 16. In the case 
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FIG. 2: Structure of the spin Hall device, a, Schematic 
image of the device. The green spheres with the red and blue 
arrows correspond to spin-up and spin-down electrons, re- 
spectively. Jc and Js are the charge and spin current respec- 
tively. The thinner red and blue arrows indicate the spin-up 
and spin-down electron motion, respectively, b, SEM image 
of the device consisting of two Py wires and a NiPd middle 
wire bridged by a Cu wire. The arrow for the magnetic field 
H represents the positive direction in the case of ISHE mea- 
surement. 



of Nio.08Pdo.92, the NLSV signal is reduced to 0.25 com- 
pared to the value without the NiPd middle wire, and 
the spin diffusion length is estimated to be 10 nm. The 
absorbed spin-up and down electrons with the orienta- 
tion sc parallel to the hard direction of the Py wires 
are scattered by Ni impurities in Pd, resulting in the 
Hall voltage Vishe- From the field variation, the overall 
change of ISHE resistance .Rishe (equal to Vishe divided 
by the charge current) can be observed. 

We show representative ISHE resistance -Rishe curves 
for Nio.08Pdo.92 as a function of the magnetic field at var- 
ious temperatures in Fig. 3c, where all curves are shifted 
for clarity. Generally -Rishe increases linearly with the 
magnetic field up to 2000 Oe and then becomes flat. This 
is consistent with the saturation of the magnetization 
of Pyl, indicated by the anisotropic magnctorcsistance 
(AMR) curve of Pyl shown in Fig. 3d. In the following 
we focus on -Rishe with saturated magnetization of Pyl 
(i.e., \H\ > 2000 Oe). For temperatures far below (10 K; 
black curve) and above (30 K; pink curve) Tc, -Rishe 
is completely flat at high enough fields, which is as ex- 
pected in the ISHE for nonmagnetic systems. The values 
of -Rishe are much larger than those for pure Pd, where 
an intrinsic mechanism should appl y 16 i 19 ' 20 and we will 
argue that this increase comes from the extrinsic effect 
of the Ni impurities. 

For temperatures near T c , in contrast, anomalous be- 
havior can be clearly seen. At 20 K just below T c 
(w 21 K), -Rishe decreases with increasing the magnetic 
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FIG. 3: Temperature dependence of the anomalous Hall and the inverse spin Hall effects of Nio.osPdo.92- a, 

AHE resistance R xy at different temperatures. The remnant magnetization M r and coercive field H c correspond to R xy at 
H = Oe and the x-intercept of R xy curve, respectively, b, Temperature dependence of M T and H c . From both plots Tc can 
be determined to be 21 K. c, Representative RISHE curves of Nio.osPdo.92 measured at 10, 20, 22 and 30 K. At 20 and 22 K, 
-Rishe at high fields is reduced and enhanced, respectively, d, A typical AMR of Pyl measured at 10 K. e, 3D plot of -Rishe 
as a function of temperature (15 — 25 K) and magnetic field, repre-sented by both height and color. 



field (\H\) up to 8000 Oe and then flattens out (see the 
red curve in Fig. 3c). Just above Tc, at 22 K, on the 
other hand, -Rishe has a positive slope versus magnetic 
field even above 2000 Oe and is not yet saturated at our 
maximum field. By making many similar scans for a 
full temperature range from 10 to 30 K, we can display 
this behavior more clearly seen in the 3D plot shown in 
Fig. 3e. The x and y axes are the magnetic field and tem- 
perature, respectively, while both the z axis and color 
represent -Rishe- When the temperature crosses Tc, a 
reduction first appears at 20 K as a light blue peak at 
negative fields or a light yellow dip at positive fields (see 
red triangles in Fig. 3c); then an enhancement appears 
at 22 K as a dark blue dip at negative fields or a red 
peak at positive fields (see blue stars in Fig. 3c). Such a 
behavior of -Rishe is observed below and above but only 
in the vicinity of T3. As already mentioned, -Rishe for 
a nonmagnetic metal should be completely flat at any 
temperature. It is thus unambiguous that the observed 
change in -Rishe of the NiPd alloys near Tc is correlated 
to the magnetic phase transi-tion of the NiPd wire (see 
Fig. 3b). On the other hand, such an anomaly near Tc 
is not observed in the temperature dependence of the 
spin absorption rate or the spin diffusion length of NiPd. 
This is because the spin diffusion length depends simply 
on the resistivity and the SO interaction of the NiPd wire 
while the ISHE signal is strongly affected by the type of 
scattering mechanisms, i.e. skew scattering, side-jump or 
band structure. 

Anomaly in ISHE of NiPd alloys near Curie tem- 
peratures. 

In order to better quantify the contribution to the 



ISHE intrinsic to the NiPd wire, we define the change in 
-Rishe as A-Rishe = [-Rishe (-ffmax) — -Rishe (—H max )]/2, 
as shown in Fig. 3c. As a maximum field H max we chose 
8800 Oe in this work. Certainly .Rishe at 22 K is not 
fully saturated even at 8800 Oe, which is different from 
-Rishe at 20 K. This is related to the fact that it is much 
more difficult to align the localized magnetic moments of 
Ni impurities above Tc compared to those below Tc even 
if high magnetic fields are applied. In Fig. 4a, we plot 
A-Rishe of Nio.08Pdo.92 (shown by red dots) as a func- 
tion of temperature. The anomaly of A-Rishe near Tc 
is clearly seen. Below 18 K, A-Rishe is almost constant 
with temperature. Only in the vicinity of Tc, however, 
it shows a dip and peak structure at 20 and 21.5 K, re- 
spectively. Since A-Rishe varies little with temperature 
except near Tc, it can be separated into the temperature- 
independent background A_Rj SHE and the anomalous 
part (5Ai?isHE, as A.Rishe = Ai?° SHE + SAR ISHE . The 
results for NiPd with 7% and 9% of Ni are plotted with 
black squares and blue triangles, respectively. The in- 
crease of Ai?j SHE with the doping concentration suggests 
the extrinsic contribution induced by the Ni impurities. 
Thus A-Rj SHE comes from both the extrinsic SHE of Ni 
impurities and the intrinsic SHE of PcU&. Despite the 
fact that A-Rj SHE depends strongly on the Ni concentra- 
tions, 5A.Rishe does so primarily via a shift to be near 
the corresponding Tc. We emphasize that such critical 
behavior was reproducible for each of at least three sam- 
ples at each Ni concentration; thus our arguments are 
based on more than 10 complete data sets, of which only 
one is shown in Fig. 3e. 

For the three concentrations, JA-Rishe are plotted as 
a function of the reduced temperature (T — Tc)/Tc in 
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FIG. 4: The anomaly in the inverse spin Hall effect 
near Tc. a, Inverse spin Hall resistance A_Rishe of Ni^Pdi-a; 
middle wires with Ni concentration x — 0.07 (black square), 
0.08 (red dot) and 0.09 (blue triangle) as a function of tem- 
perature. The anomaly of ATIishe for each Ni concentra- 
tion appears only near its Tc- b, The anomalous compo- 
nent <5Ai?isHE for Ni concentration equal to x — 0.07 (black 
square), 0.08 (red dot) and 0.09 (blue triangle) plotted as a 
function of normalized temperature (T — Tc)/Tc- For com- 
parison we also show the uniform second order nonlinear sus- 
ceptibility X2 nl (black solid line). 



Fig. 4b. Note that each curve corresponds to one full 
data set: for example the red line represents a reduction 
of all of Fig. 3e. Although the three devices have different 
Tc and Ai?j SHE , the anomalous behavior appears to be 
universal near Tc: the three curves almost scale onto one. 
The amplitudes of the dip and peak are about 5 fifl. 
The anomalous behavior disappears when the reduced 
temperature differs from Tc by 10%. These results are 
obviously different from the case of the AHEi£ where the 
anomaly of the Hall resistivity appears only below Tc . 

Discussion 

What is the origin of the anomaly of SARishe observed 
only in the vicinity of Tc? Resistive anomalies near mag- 
netic critical points Tc have been studied for a long time. 
The anomaly near Tc in the longitudinal resistivity ap- 
pears within the first Born approximation from quadratic 
fluctuations ((Mi - (Mi))(Mj - (My))) in the moments 
Mi and Mj of localized spins at sites i and j where con- 



tributions come only when i and j are within a certain 
cut-off distance^. The resistivity then depends on a lin- 
ear local magnetic susceptibility Xo° c which is a summa- 
tion of short-range fluctuations up to this distance. To 
simplify notation we shall denote the sum of such short 
range correlations as ^2((M — (M)) 2 ). omitting the site 
subscripts, and we shall use similarly abbreviated nota- 
tions for higher order moments. As mentioned in the 
introduction, the anomaly of pn in the AHE near Tc was 
explained by Kondoi^. Note that the quantities r\ and T2 
in ref. 13 correspond to a first order nonlinear suscepti- 
bility xi° c ~ — (-^O) 3 ) an d a second order nonlin- 
ear susceptibility X2° c 2((-^ — (^)) 4 )i respectively. We 
change notation for simplicity, and also to emphasize that 
Kondo makes a strictly local approximation, whereas we 
extend the theory to include nearby correlations between 
atomic moments. The anomaly of pn in the AHE ap- 
pears from the first order fluctuation of local moments 
Xi° c - The second order spin fluctuation X2° c is also de- 
fined in ref. 13 but made no contribution to the AHE. 
This will be explained below; we first note that the Xi° c 
and X2° c terms appear in the s-d Hamiltonian to the same 
order with respect to the spin-orbit coupling constant A 
of the localized moment. The term in Xi° c changes sign 
for spin-up or down electrons, while X2° c does not. Thus 
the scattering amplitudes for spin-up and down conduc- 
tion electrons are different and proportional to xi° c + X2° c 
and — xi° c + X2° c respectively- 1 ^. As illustrated in Fig. I, 
in the case of the AHE configuration, spin-up and down 
conduction electrons are scattered to opposite directions. 
As a result, pn is proportional to the difference of scat- 
tering amplitudes between spin-up and down electrons, 
i.e., xl° c - O n the other hand, in the 1SHE configura- 
tion, both spin-up and down electrons are scattered to 
the same side. Unlike the AHE, the obtained -Rishe is 
proportional to the sum of scattering amplitudes between 
spin-up and down electrons, i.e., X2° c - The effect of such 
terms has been hidden for 50 years 1 ^. However, it is nat- 
urally expected that the contributions related to X2° c can 
survive if the incident current is the spin current, which 
is exactly the case for the ISHE. 

To make a comparison between theory and experiment 
we should, in principle, calculate X2° C j which is integrated 
over the Fermi surface. This X2° c should be proportional 
to the experimental <5Ai?isHE- However, the integration 
depends essentially on the Fermi wave vector, which acts 
as a cut-off to any divergence. To avoid including mate- 
rial dependent parameters, we compare to a quantity, the 
uniform second order nonlinear susceptibility x™ 1 which 
is independent of cut-off. xlf" is the response of the mo- 
ments on all the sites to a uniform magnetic field, and 
is the sum of correlations over all distances without any 
cut-off (see Supplementary Discussion). In Fig. 4b we 
superimpose on the experimental data a mean field cal- 
culation of X2 nl - It diverges with opposites signs above 
and below Tc . We interpret the dip and peak observed in 
the experimental data as a smearing of such a divergence. 
In a full calculation the theoretical spin Hall resistance 
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should also be smeared: the divergence is cut-off, just as 
for the longitudinal resistance^!. To confirm the above 
scenario, we have performed Monte Carlo simulations for 
X 2 oc where the contribution of nearest-neighbor correla- 
tions is taken into account. As can be seen in Supplemen- 
tary Figure SI, the smearing of the divergence of X2° c at 
Tq is well-reproduced. Thus we argue that neither the 
full theory, nor the experiments should diverge at To 

Finally we note that the anomaly of pn in the AHE in 
pure Feii and N i 11 ' 12 may alternatively be explained by 
Karplus and Luttinger's theory which is based on band 
effects 2 ^. In this theory pn would be proportional to 
p xx (M) where p xx is the longitudinal resistivity and (M) 
is the averaged magnetization. Since (M) is zero above 
Tc , the sign change of i?isHE below and above Tc seems 
to preclude explanation by such band effects. 

In conclusion, we have demonstrated a characteristic 
behavior of the spin Hall effect at the magnetic phase 
transition of the weak ferromagnetic NiPd near its Curie 
temperature To It can be explained by the fluctuation 
contributions to skew scattering via spin-orbit interac- 
tions, and appears as the second order nonlinear spin 
fluctuation X2° c - Considering the small volume of NiPd 
wire (10 -16 cm 3 ) and the magnetic moment of NiPd, 
the total magnetic moment involved in the experiment 
is extremely small: less than 10~ 14 emu. Such a tiny 
moment is far below the sensitivity of magnetometer (for 
instance, 10~ 8 emu for SQUID). The successful extrac- 
tions of the contribution of higher order spin fluctua- 
tions in our experiments prove the extremely high sen- 
sitivity of our method. This technique could be very 
useful to understand more complex spin textures such 
as spin-glasses and spin-ice. For those systems, the lin- 
ear magnetic susceptibility Xo" 1 has been well-studied so 
far but there have been almost no experimental studies 
of the higher order spin fluctuations where some impor- 
tant physics might be still hidden. The present method 
constitutes a novel approach to detect high order nonlin- 



ear fluctuations with an extremely high sensitivity, which 
is impossible to achieve by a direct measurement. This 
will broaden the application of the SHE from potential 
devices to a means of measuring higher order spin fluc- 
tuations, particularly in the critical region. 

Method 

Sample preparation. 

The spin Hall devices were fabricated on a thermally 
oxidized silicon substrate using electron beam lithogra- 
phy on poly-methyl-methacrylate (PMMA) resist and a 
subsequent lift-off process. Here we used a lateral spin 
valve structure which consists of two Py wires (30 nm 
thick and 100 nm wide) with a separation distance of 
1 pm and a NiPd middle wire (20 nm thick and 100 nm 
wide) bridged by a Cu wire (100 nm thick and 100 nm 
wide) as is detailed in Fig. 2 in the main text. A Pd wire 
was first deposited by electron beam evaporation using a 
99.999% purity source. In order to form the NiPd alloy 
wire with a Ni content of 7%, 8% and 9%, we implanted a 
suitable amount of Ni + into the Pd wire using an Ion Im- 
planter (NH-20SR) with an acceleration voltage of 30 kV. 
As is mentioned in the main text, the prepared NiPd al- 
loy is not perfectly homogenous and nickel is distributed 
over ~ 10 nm below the Cu/NiPd interface in 20 nm 
thick Pd wires. This nickel distribution would lower the 
Tc compared to bulk NiPd alloys. However, the spin dif- 
fusion length of NiPd of about 10 nm corresponds to the 
maximum distribution of Ni in Pd. That is why we could 
still observe the spin fluctuation of Ni moment. The two 
Py wires were deposited by electron beam evaporation, 
while the Cu bridge was fabricated by a Joule heating 
evaporator. Prior to the Cu evaporation, a careful Ar ion 
beam milling with a beam current of 12 mA and beam 
voltage of 600 V was performed for 30 seconds to clean 
the surfaces of Py and NiPd wires and to obtain highly 
transparent contacts. 
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